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Abstract 

We compute one-loop effective action of magnetic monopole in three-dimensional 
electrodynamics of massless bosons and fermions and find that it contains an in- 
frared logarithm. So, when the number of massless matter species is sufficiently 
large, monopoles are suppressed and in the weak coupling limit charged particles are 
unconfined. This result provides some support to gauge theories of high-temperature 
superconductors. It also provides a mechanism by which interlayer tunneling of exci- 
tations with one unit of the ordinary electric charge can be suppressed while that of 
a doubly charged object is allowed. 
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Gauge theories of high-temperature superconductors assume that in certain 
planar electronic systems spin and charge are separated and the resulting new quasi- 
particles interact via abelian gauge forces. This idea is seemingly in contradiction 
with three-dimensional confinement due to magnetic monopoles @. One may try 
to resolve the contradiction by assuming that some of the charged excitations are 
gapless. Then it is possible that quantum effects of these excitations are sufficiently 
strong in the infrared to modify the interaction between monopoles in such a way that 
confinement of charges is lost. To see how this can happen, consider the relativis- 
tic version of the problem - three-dimensional quantum electrodynamics of massless 
bosons or fermions. A simple calculation shows that one-loop contribution of massless 
charged particles to the gauge field propagator causes its small momentum behavior 
to change from the usual 1/p 2 to l/[p|. The bilinear part of the corresponding effec- 
tive action for the gauge field calculated on the monopole configuration then produces 
a logarithm of the total size of the system. Indeed, the field strength of monopole is 
Fij ~ €ijkXk/r 3 , therefore, 

/ d * xF «iw Fij ~ ' (1) 

where R and a are infrared and ultraviolet cutoffs, respectively. Essentially the same 
calculation appears in ref. ||. The coefficient of the logarithm increases proportionally 
to the number N of massless species. This suggests that in the presence of massless 
charged particles, at least when N is large enough, monopoles are suppressed via a 
version of the "infrared catastrophe" and charged particles are unconfined. 

The reason why eq. ([I]) is not sufficient to determine the fate of magnetic monopoles 
in the presence of massless charged particles even in the large N limit, is that the 
interaction of a charged particle with monopole has no small parameter, so the bi- 
linear part (0) of the effective action is in no way distinguished relative to terms 
containing more powers of the gauge field. To make a reliable conclusion, we need 
all these terms. On the other hand, because the gauge propagator in the large N 
limit is of order 1/N, the large N limit suppresses higher-loop contributions to the 
effective action. Therefore, in this limit the problem reduces to calculation of the 
one-loop determinants of massless bosons and fermions in the monopole background. 
For this calculation, we chose to proceed with relativistic three-dimensional particles. 
Precise dispersion laws for quasiparticles that may occur in real electronic systems are 
unknown at present. The result and the main steps of our calculation are presented 
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below. The result confirms the presence of an infrared logarithm in monopole's effec- 
tive action both in bosonic and fermionic cases. We thus show that single monopoles 
are suppressed by an "infrared catastrophe" in the presence of a sufficient number 
of massless matter fields. We cannot state at present what exactly this "sufficient" 
number is because the answer to this question lies outside the region of validity of 
the large N approximation. 

Consider now the weak coupling limit of three-dimensional electrodynamics when 
the dimensionful gauge coupling e 2 is much smaller than the ultraviolet scale M ~ a -1 
at which internal structure of monopole becomes essential. Then, in the absence of 
massless matter, monopoles and anti-monopoles would form a dilute gas. When the 
number of massless matter species is sufficiently large, the infrared logarithm causes 
monopoles and anti-monopoles to assemble into "molecules" - pairs of typical size 
d that is much smaller than the average distance between the pairs. These pairs 
interact by a short range potential of order (d 2 /r 2 )ln(d/r), so it is natural to expect 
that charged particles are unconfined. This picture provides some support to the idea 
of new gauge interactions in planar electronic systems. At the end of this paper we 
discuss some further applications of our results. 

The one-loop contributions to monopole's effective action from a single charged 
bosonic field and a single charged fermionic field are respectively of the form 

S B 1] = Trln(— D 2 ) - Trln(-<9 2 ) , S ( F ] = -Trln(aD) + Trln(ad) , (2) 

where D are covariant derivatives and a are Pauli matrices. Eq.([|) requires both 
ultraviolet and infrared regularizations. We compute not expressions (@) directly but 
rather 

S$\R) = TrlmM b - TrlmM bo , S^(R) = ~ (Trln(-M 2 F ) - Tr\n(-M 2 F0 )) , 

(3) 

where 

Mb = -^(r 2 + R 2 )D 2 (r 2 + R 2 ) , -M 2 F = -^(r 2 +R 2 )(aD) 2 (r 2 + R 2 ) , (4) 

and Aiso and — Ai 2 F0 are obtained analogously from the operators — d 2 and —{ad) 2 
of the vacuum sector. This replacement is similar to the one used by 't Hooft in his 
four-dimensional instanton calculation 0]. If the effective action were not infrared 
sensitive, the additional factors of (r 2 + R 2 )/R 2 would cancel between vacuum and 



3 



non- vacuum contributions in (|3]). In the present case, we intend to show that the 
effective action is infrared sensitive. In this case eq.(|3]) provides an infrared regu- 
larization of eq.(|2]), R being the regulator radius. In what follows we measure all 
distances in units of R, hence, we set R — 1. 

The eigenvalue equations for operators (§) are 

= , (5) 




, ,>,- , - ■ 0, (6) 

where A stand generically for the eigenvalues. In bosonic eq.(|5|), radial and angu- 
lar variables are separated by ^/b = V ; ( r )^g,z,m(^, 0), where Y q ^^ m are the monopole 
harmonics of ref. ||. One gets the radial equation 
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a(a + l) 4A 
+ 



^•2 f' 1 j .y*2\2 



i> = o , (7) 



where a — [(I + 1/2) 2 - q 2 ] 1/2 - 1/2, / = |g|, |g| + 1, and the multiplicity of the 
eigenvalue A is (21 + 1). Parameter q assumes integer and half-integer values as a 
consequence of the Dirac quantization condition. All our results depend only on \q\, 
so in what follows we take q > 0. In fermionic eq. (]?]), the variables are separated 
by using either of the three angular dependences Q^l, f] m introduced in ref.||. 
We find that in all three cases the resulting radial equations have the same form as 
eq.(0) but with different values of a. For angular dependence ^\ a = jj — 1, where 
A* = [(j + 1/2) 2 — q 2 } 1 / 2 and j = q + 1/2, q + 3/2, while for angular dependences 
and 7], a = fx with j = q + 1/2, q + 3/2, ... for £( 2 ) and j = q — 1/2 for 77. In all cases 
the multiplicity is (2j + 1). In three dimensions, fermionic wave functions (or, more 
precisely, wave sections 0, |6| ) in © are doublets. This doublet structure is carried 
by the angular dependences, so both for bosons and fermions radial functions ip in 
eq.(|7]) are one-component objects. Therefore, unlike the scattering problem in (3+1) 
dimensions 0, our calculation does not require any special treatment of j = q — 1/2 
fermionic modes. 

We can treat bosonic and fermionic cases simultaneously using eq.(|7|) if we adopt 
the following notation 

a = [(j + 1/2) 2 - q 2 } 1 ' 2 - k - 1/2 , j = q + K , q + k + 1, ... . (8) 
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Then, k = corresponds to bosons, k = 1/2 to fermions with angular dependence 
a nd k = — 1/2 combines fermions with angular dependencies ^ 2 ^ and 77. Results 
for the vacuum sector are obtained by substituting q = 0. Though such substitution 
leads to an unphysical value j = —1/2 for k = —1/2, eigenvalues corresponding to 
this unphysical value do not participate in the fermionic trace in (§) because of the 
vanishing multiplicity factor (2j + 1). 
By the change of variables 

${r) =r a (l + r 2 )- Q ~ 1/2 0(x) , x=(l+r 2 )" 1 , (9) 

eq.(|7|) is converted into a hypergeometric equation. The resulting eigenvalues are 

A n = (n + a + l/2)(n + a + 3/2) , n = 0, 1,... . (10) 

We still need an ultraviolet regularization for the traces in eq.(||). A convenient 
one is provided, again in parallel with 't Hooft's calculation [[|, by two Pauli-Villars 
regulators with masses Mj and metrics satisfying J2i e i — ~ 1> J2i e«M 2 = 0, 
i = 1,2. Regularized traces that we will need are 

00 00 

TrlnM(M i; n)= £ (2j + 1) £ £ eJn[( S + a) 2 + /i 2 ] , (11) 

j=g+K s=l 1=0,1,2 

where we defined e = 1, M = 0, and /i 2 = M 2 — 1/4, 2 = 0, 1, 2. The effective action 



in the fermionic case is obtained from the half-sum of traces (11) with k = ±1/2 which 
are in fact related to each other in a simple way. 

Because the effective action of monopole is dimensionless, it can depend only on 
products of infrared and ultraviolet regulator parameters, MiR. In the system of units 
where R — 1, we are then interested in the dependence of the effective action on Mj 
in the limit when Mi are large. Non- regulator terms in eq.flTTD cannot produce such 
dependence. Applying the Euler-Maclaurin formula to regulator terms in eq. (O) we 
get 

T \n[(s + a) 2 + M 2 } = [ A ds\n[(s + a) 2 + M 2 } 

s=l J ° 



+ (hn[(s + a) 2 + M 2 } + l- ■ s + " 



where A 3> M. The remainder in eq. (|12]) gives a convergent contribution of order 1/M 
when summed over j with the multiplicity factor (2j + l) and hence may be neglected. 
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Terms divergent at A — > oo as well as those proportional to M 2 get cancelled when 
all regulator and non-regulator contributions are added together and we obtain 



TrlnA4(Mj; k) = £ (2j + 1) 



x 



j=q+K 

E e i (-(« + l/2)ln(« 2 + ^) + 2^ ; (^-arctan— ] -- ° 2 ) +C(J) 

" (13) 

where C(j) denotes terms independent of Mj. 

The dependence of eq. (|13"1) on Mj can be found by the following method. We 
decompose each sum over j into two - one running from q + k, to some value J — 1 
such that q J Mi, another running from J to a cutoff A' 3> M. As in the 
sum over s before, dependence on the cutoff will disappear when all regulator and 
non-regulator terms are added together. The number J is integer or half-integer when 
q + k is integer or half-integer, respectively. Now, in the region q + k < j < J — 1 
we can neglect j compared to Mj while in the region J < j < A' we can use the 
Euler-Maclaurin formula. At J < j < A' we can also use the expansion 

a = ai - q 2 (2j + l)- 1 + 0(r 3 ) , a x =j-K. (14) 

For example, 

A' 

£ (2j + l)(a + l/2)ln(a 2 + M 2 ) 

j=q+K 

J-1 A' 

= ]T (2j + l)(a + l/2)lnM 2 + £(2j + l)( ai + l/2)ln(a 2 + M 2 ) 

j=g+« j=J 

- g 2 f [h(o? + M 2 ) + a ^+ 2 1} ) + 0(M- 2 ) + 0( J- 1 ) . (15) 

It turns out that other terms in eq . C[3~3|) do not produce contributions in either bo- 
son or fermion determinant that distinguish between monopole and vacuum sectors. 



Proceeding with eq . (|T5|) , we finally obtain the one- loop effective action of a monopole 
with monopole number q in the presence of massless bosons and fermions up to terms 
independent of R: 

S { ^ F (R) = K B , F (q)\nM 2 R 2 + 0(R°) , 



K B (q) = Jim 

J— »oo 



J2(2j + l)[(j + 1/2) 2 - q 2 f 2 - + + q 2 J 

3 6 



(16) 
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Q 






1/2 


0.0968 


0.0151 


1 


0.2266 


0.1730 


3/2 


0.3850 


0.4358 


2 


0.5682 


0.7852 



Table 1: Coefficients of \nR 2 in the one- loop effective action of a monopole with 
monopole number q in the presence of massless bosons and fermions. 



K F (q) 



lim 

J— >oo 



£ (2 J + l)[( J + l/2) 2 -,r! 

j=q+l/2 



211/2 



2 t3 1 t 
-J + - J 

3 6 



q 2 J 



Q 
2 



(17) 



We assumed that both regulator masses are of the same order, Mi ~ M. The limits in 
eqs.(|T6|)-(|TT|) were done by computer. The results for a few values of q are presented 
in Table 1. 

When there are N species of particles of given type, the corresponding numbers 
from Table 1 should be multiplied by N. For sufficiently large N, logarithms of 
R coming from boson or fermion determinants will overpower 31nMi? that comes 
with the opposite sign from the volume factor. Thus, when there is a large number of 
massless matter species, monopoles are suppressed. For small N, we cannot draw any 
conclusions from the present work because the one-loop calculation is not a reliable 
guide in this case. Possibly, numerical simulations can help to solve the problem for 
small N. 

Let us now state some applications of our results. They show that non-confining 
abelian gauge interactions are possible in (2+1) dimensions when gapless excitations 
are present, thus providing some support to the idea of new gauge interactions in 
planar electronic systems. In a somewhat different interpretation, our results provide 
a mechanism by which interlayer tunneling of excitations with one unit of the ordi- 
nary electric charge is suppressed while that of a doubly charged object is allowed. 
Recent work J7| shows that a viable theory of high-temperature superconductors can 
be constructed if such mechanism is assumed to exist. Let us postulate that a planar 
system describing one layer in a layered material supports quasiparticles that carry 
magnetic flux with respect to the new gauge field. There are two varieties of such 
quasiparticles corresponding to positive and negative fluxes, respectively. Assume 
further that these quasiparticles carry ordinary electric charge (" holons" ) which has 
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the same sign both for positive and negative fluxes. Events of interlayer tunneling 
of these quasiparticles are described by monopoles and anti-monopoles in three di- 
mensions. If there are also excitations of another sort ("spinons") that carry charge, 
rather than flux, with respect to the new gauge field and are gapless, tunneling of a 
single flux can be suppressed by their infrared effects as discussed above, while a pair 
of positive and negative fluxes can tunnel freely. 
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Rubakov for discussions on fermion-monopole interactions. It is a pleasure to acknowl- 
edge the hospitality of the Aspen Center for Physics where this work was started. The 
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